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Voltage-dependent ion channels determine the electric properties of axonal cell 
membranes. They not only allow the passage of ions through the cell membrane 
but also contribute to an additional charging of the cell membrane resulting in the 
so-called capacitance loading. The switching of the channel gates between an open 
and a closed configuration is intrinsically related to the movement of gating charge 
within the cell membrane. At the beginning of an action potential the transient 
gating current is opposite to the direction of the current of sodium ions through the 
membrane. Therefore, the excitability is expected to become reduced due to the 
influence of a gating current. Our stochastic Hodgkin-Huxley like modeling takes 
into account both the channel noise - i.e. the fluctuations of the number of open 
ion channels - and the capacitance fluctuations that result from the dynamics of the 
gating charge. We investigate the spiking dynamics of membrane patches of variable 
size and analyze the statistics of the spontaneous spiking. As a main result, we find 
that the gating currents yield a drastic reduction of the spontaneous spiking rate for 
sufficiently large ion channel clusters. Consequently, this demonstrates a prominent 
mechanism for channel noise reduction. 

PACS numbers: 05.40.-a, 87.16.-b, 87.19.Nn 



I. INTRODUCTION 

Following the study of Hodgkin and Huxley most of the models of axons have treated 
the generation and propagation of action potentials using deterministic differential equations. 
Since the work of Lecar and Nossal it has become increasingly evident, however, that 
not only the synaptic noise but also the randomness of the ion channel gating itself may 
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cause threshold fluctuations in neurons 0, 4|. Therefore, channel noise which originates in 
the stochastic nature of the ion channel dynamics should be taken into account [3|, |j] . For 
example, in mammalian ganglion cells both the synaptic noise and the channel noise might 
equally contribute to the neuronal spikes variability [5]. Due to finite size, the origin of 
the channel noise is basically due to fiuctuations of the mean number of open ion channels 
around the corresponding mean values. Therefore, the strength of the channel noise is 
mainly determined by the number of ion channels participating in the generation of action 



potentials. Channel noise impacts, for 



and the spiking rate itself 
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e, such features as the threshold to spiking 



121], the anomalous noise-assisted enhancement 



of transduction of external signals, i.e. the phenomenon of Stochastic Resonance 
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,122|,123|,|24| and the efficiency for synchronization [251]. Interestingly 
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enough, there exist optimal patch sizes for which the spike production becomes more regular. 
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23l . |26| . The objective of this 



or the response to external stimuli optimizes 1^, 
work is to investigate the spontaneous spiking of cell membrane patches when gating charge 
effects are also considered. 

When an ion channel opens or closes, an effective gating charge is moved across the 
membrane [27 1 . This motion creates the so-called gating current which is experimentally 
measurable |28l. l29l. l30l]. The influence of gating currents was not explicitly considered in the 
orig inal Hodgkin-Huxley (HH) model [l]. In 1975 Hodgkin 31] and, independently, Adrian 



321 ] , flrstly inquired theoretically into the influence of the ion channel density on the velocity 
of the action potential propagation along the squid giant axon by taking into consideration 



the gating currents o: 
ion channels, see also 



sodium ion channels (via an effective capacitance loading caused by 



331]). Remar 
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the signal velocity is maximal 
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they found an optimal ion channel density for which 



35, 



361 ]. This led to the emergence of the Hodgkin's 



maximum velocity hypothesis, i.e. that the living species can adopt optimal densities of ion 
channels to maximize the rate of neuronal information transduction. 

In this work we investigate the influence of both the gating currents and the channel 
noise on the spontaneous spiking. At the initiation of a spike, there occur two competitive 
currents in the presence of ion channels: First, an ionic current, which is caused by the 
influx of positively charged sodium ions through open sodium ion channels into the cell, and 
second, a gating current that flows in the opposite direction during the transition of sodium 
ion channels from closed to open state. The latter one is due to the movement of the positive 
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gating charge from the inside of the cell membrane to the outside. This current actually 
precedes the first one, and terminates when the channel is in its open state. However, due 
to spontaneous fluctuations between the channel's states, a stochastic component of gating 
current also emerges. Moreover, there are also gating currents due to p otassium channels 



which have been neglected within the macroscopic, noise- free model |31l. I32l. I34J. I35[|. These 
polarization currents of membrane-bound charges can substantially contribute to the overall 
polarization current across the membrane, which can be considered as an electrical capacitor. 
Therefore, the influence of gating currents on the excitability of small ion channel clusters 
is a priori not clear and necessitates further research. 

In section [III we show how channel noise can be taken into account within a Langevin 
equation approach. The next section IIIII deals with the extensions to the standard HH 
model with respect to gating charge effects. We discuss the influence of gating charge on 
the dynamics with respect to excitability due to a constant stimulus and as well as on the 
deterministic spiking period. Thenceforth, in section [TV] we consider the situation where 
both channel noise and gating charge are considered and demonstrate the influence on the 
mean interspike interval. The effect of intrinsic coherence resonance in the presence of gating 
currents is discussed in section |Vl 

II. A STOCHASTIC GENERALIZATION OF THE HODGKIN-HUXLEY 

MODEL 

According to the standard HH model the dynamics of the membrane potential V follows 

C^V + GK{n) {V - E^) + GNalm, h) {V - E^,) + Gi^ {V - El) = /ext(t) , (1) 

with the potassium and sodium conductances given by: 

G'KH=^7^^n^ G^,im,h)=g^:-m'h. (2) 

In Eq. ([1]), V denotes the membrane potential measured throughout this work in [mV] 
and G = 1/iF/cm^ is the capacity of the cell membrane. The time t is scaled in ms. 
Furthermore, -ENa = 50 mV, En = — 77mV and -El = — 54.4mV are the reversal potentials 
for the potassium, sodium and leakage currents, correspondingly. The leakage conductance is 
assumed to be constant. Gl = 0.3mS/cm^, and g^^^ = 36mS/cm^ and g^^^ = 120mS/cm^ 
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denote the maximal potassium and sodium conductances, when all ion channels are open. 
Iext{t) indicates an external current stimulus. Eq. ([1]) is nothing but a Kirchhoff law for an 
electrical circuit made of membrane capacitor and variable nonlinear conductances assuming 
that the conductances of open ion channels are Ohmic, i.e. all the nonlinearity comes from 
their gating behavior. 

The gating variables n, m and h, cf. Eq. ([1]) and Eq. ([2]), describe the mean ratios of 
the open gates of the specific ion channels. Assuming gate independence, the factors ra^ 
and h are the mean portions of the open ion channels within a membrane patch. The 
dynamics of the gating variables are determined by voltage-dependent opening and closing 
rates a^{V) and (3^{V) {x = m, h, n) taken at T = 6.3 °C. They read Q, Q, 42]: 



1^ + 40 

"-^'^^ = °- S-exp{-(V + 40)/10} ' 

(3UV) = 4 exp {- {V + 65)/18} , (3b) 
ahiV) = 0.07 exp{-{V + 65)/20} , (3c) 

^^^'^^= l + exp{-(V + 35)/10} ' ^'^^ 

V + 55 

aJV) = 0.01 — -— ^ , (3e) 

^ ' 1 -exp {- (1^ + 55) /lO} ' ^ ^ 

/3„(y) = 0.125 exp {-(V + 65)/80} . (3f) 

The dynamics of the mean fractions of open gates reduces in the standard HH model to 
relaxation dynamics: 

^x = a^{V) {1 - x) - PxiV) X, x = m,h,n. (4) 

Such an approximation is valid for very large numbers of ion channels and whenever fluctu- 
ations around their mean values are negligible. 

However, each channel defines in fact a bistable stochastic element which fluctuates be- 
tween its closed and open states. The same is valid for the gates which assumed to be 
independent in the HH model. The number of open gates undergoes a birth-and-death-like 
process. The corresponding master equation can readily be written down |l^ . The use of a 
Kramers-Moyal expansion in that equation results in a corresponding Fokker-Planck equa- 
tion which provides a diffusional approximation to the discrete dynamics. The corresponding 
Langevin equations then reads: 

^x = a^iV) {1 - x) - I3.j:{V) X + ^a:it), x = m,h,n. (5) 
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Here, Cx{t) are independent Gaussian white noise sources of vanishing mean. For an excitable 
membrane patch with A^^a sodium and A^k potassium ion channels the noise correlations 
assume the following form: 



(6a) 
(6b) 
(6c) 



{Ut)Ut')) = ]vL K(V^) (1 - m) + PUV) m] Sit - t') , 
(UW)) = jt i^^) (1 - ^) + f^hiV) h] 6it - t') , 
{Ut)Ut')) = i4 K(V) (1 - n) + Pr^iV) n] S{t - t') . 

The stochastic Eq. ([5]) replaces Eq. (jl]). Note that the correlations of the stochastic forces 
in these Langevin equations contain the corresponding state-dependent variables and thus 
should be Ito-interpreted because it is derived in a diffusional approximation of a jump 



process which is intrinsically a white noise problem 
white noise limit of a colored noise problem 40l |. 
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38, 



391], and does not result from a 



Within the approximation of homogeneous ion channel densities, pNa = 60 yum~^ and 
Pk = 18 /im~^, the ion channel numbers are given by: A^Na = PNa'S', A^k = PkS, with S 
being the size of the membrane patch. The number of ion channels, or the size of the excitable 
membrane patch S, respectively, determines the strength of the fluctuations and thus the 
channel noise level. With decreasing patch size, i.e. decreasing number of ion channels, the 
noise level caused by fluctuations of the number of open ion channels increases, cf. Eq. (l6a|l - 
fl6cl) . Within this modelling the rate for spontaneous spiking depends monotonously on the 



patch size 



21 



25|. 



III. GATING CHARGES AND CURRENTS IN A MODIFIED 
HODGKIN-HUXLEY MODEL 



Our starting point for the derivation of the gating charge values are the opening and 
closing rates, Eq. ( |3al) - (l3fj) . of the standard HH model. Conformational changes of the 



ion channel in order to switch between open and closed configurations come along with 
movement of gating charge jl, Assuming Arrhenius-like dependences, the transition 



rates are 



3: 



a{V) = aoexp 
(3{V) = /3oexp 



qz ■ V 
-q{l-z)-V 



(7a) 
(7b) 



6 



containing the voltage-independent parts ao, /3q, the gating charge q, and the asymmetry 
parameter z with z G (0, 1). The two latter parameters can be deduced by comparison of 
Eqs. (I7aD-(l7bl) with Eqs. (E^h 



large negative voltage values [41 



3f|) in the high activation barrier limit, being assumed at 
42|], yielding 



aniV) oc exp {+ V/10} , f3niV) oc exp {- V/80} , (8a) 
ahiV) (X exp {- V/20} , P^iV) cx exp {+ V/IO} , (8b) 
aUV) cx exp {+ V/10} , PUV) oc exp {- V/18} . (8c) 



The corresponding gating charges then follow as: 



qn = ksT 1000 ■ + ^) = 2.709e , (9a) 
qh = ksT 1000 ■ (-2^ + To) = -3.612e , (9b) 
q^ = ksT 1000 ■ (^ + ^) = 3.746e , (9c) 

ar, q^ and g„ correspond to 



271] from the interior side of 



with the elementary charge e = 1.6022 • 10^ A ■ s. In particu 
motion of effective positive charges during the gate openings 
cell membrane to the exterior one (outward motion). The total gating charges are 

gn-charge = ^ 4 TVk , (10a) 
gm-charge = x 3 iV^a , (10b) 
^h-chargc ^^^xliVNa. (lOc) 

To obtain the gating currents one has to multiply the gating charges with time-derivatives 
of the corresponding gating variables. Adding the corresponding current densities (per unit 
area) to the equation for voltage variable ([1]) we find 

/ext(t) = + g^l^nv" h{V- \/Na) + gT^'n^ iY - ^k) + 9i{V - Vi) 



dt 

+PNa3^g„^ + PNaf g/. + PK4f g., (11) 



where px are the specific ion channel densities. The dynamics of the gating variables are 
given by Eq. ([5]) and Eq. (!6al) -( l6cll . Taken together, all these equations constitute a stochastic 
generalization of the HH model which accounts for gating current effects being considered 
in the present work for the first time. 
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FIG. 1: The bifurcation diagrams for the original deterministic HH model in ([T]) (black lines) 
and the deterministic generalized model (grey lines) in (llip which takes gating charge effects 
into account, are plotted versus the external constant driving current /ext- The continuous lines 
correspond to stable solutions, dotted lines refer to the unstable solutions. In Fig. (a) the minimal 
and maximal membrane potentials are plotted. In Fig. (b) the period of the spiking is depicted. 
A horizontal line at zero indicates the existence of a non-spiking solution. 
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Deterministic case: neglecting channel noise 

In order to study the influence of gating cfiarge on tlie generation of action potentials 
in more detail we consider the case of constant current driving, i.e. Icxt{t) = const. Fig. [1] 
depicts the bifurcation scenario for spike occurrence in the original HH model and in the 
generalized model which takes gating charge into account. As, for the moment, we want 
to neglect channel noise (which corresponds to the limit of infinite patch size or infinite 
numbers of ion channels), the gating dynamics, Eq. [5l reduces to that in Eq. (j4]). 

While the existence of the fixed point solution is not affected by the gating charge effects, 
both its stability range and the emergence of the oscillatory spiking solution are shifted 
towards larger values of the external driving current, cf. Fig. [1] (a). This again verifies the 
inhibitory functionality of the gating charge. For the present case, the critical current value 
at which the fixed point loses stability is Ip^t ~ 10.81 /lA/cm^, whereas in the standard HH 



model it occurs at /gxt ~ 9.76yuA/cm^ 43|, |4J]. Yet, before this happens in both models 



in ([T]) and in (fTTl) . a periodic, spiking solution emerges via a subcritical Hopf bifurcation. 
Then both the fixed point solution and the limit cycle, i.e. the spiking solution, coexist 
for a certain range of driving currents. The periods of the solutions are depicted in Fig. [1] 
(b). The spiking period vs. the driving current is also shifted towards higher values of the 
driving current. 

Fig. [1] demonstrates that the gating charge does not affect the rest potential. It impedes, 
however, the system's excitability. For a given Jgxt the threshold for excitation is increased, 
meaning that the excitability in the case of a gating charge dynamics becomes decreased 
as compared to the case when no gating charge effects are present. In addition, the time 
period of an oscillatory spiking solution is increased. 



Deterministic case: increment of the membrane time constant 

At the beginning of an excitation the m gates open, therefore ^ > 0. Consequently, the 
m-gating current is positive (corresponding to an outward motion of positive charges), i.e. 
-^NaT^^*™^ = '^Nj^aQm Iff > {Qiu > 0). The opening of sodium channels enables transport of 
positively charged sodium ions into the cytoplasm (implying a negative current by standard 
convention). Thus, these two currents (the gating or the polarization current stemming 
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from the ion channels and the ion current) are oppositely directed. At the very beginning 
of the spike initiation (the gating variable m is ca. 0.1 at the threshold) the absolute 
value of the gating current is around 10 % of the ionic current. The reason is that the 
gating current is proportional to the m variable, whereas the averaged ion current scales 
with m^. Moreover, the gating current enhances the overall polarization current across 
the membrane. As a result, the membrane time constant becomes effectively enhanced 
as it can be deduced from the following reasoning: To study the spike initiation, only 
the fastest gating variable m(t), which rapidly adjusts to the steady state value, m^iV) = 
Oimiy)/ [amiy) + Pmiy)] should be taken into account. Note also that n{t) and h{t) exhibit 
a slow dynamics and their values can be taken at the rest potential V^est = —65.0 mV. It then 
follows, using ^ 'i'^^^v) ^i^Q influence of gating currents can be accounted 

for by a voltage-dependent contribution to the membrane capacitance. Consequently, taking 
this m-gating charge contribution into account, the effective membrane capacitance reads 
C^g{V) =C + Cg{V), where 

Cg{V) = 3pj^^qm — — = 3pNagm"^oo(V^) [l-moo(V^)J ^ m ^ (12) 



This latter relation assumes a simple and insightful form in the approximation of Eqs. flTal) . 

CgiV) ^ 3pNa7%moo(V^)[l - m^iV)]. (13) 

This expression evidences that Cg{V) is maximal when the sodium m-gates are half-open, 
i.e. moo = 1/2. For qm given in Eq. ( l9c|) and with T = 6.3 °C, we have Cg ~ 2.8 • 
10~^^ PNa"^oo(^)[l — '^oo(^)] F. Furthermore, for pNa = 60 /xm~^, and for nioo ~ 0.05 at the 
resting potential, Cg(Vrest) ~ 0.085 /iF/cm^, i.e. the additional capacitance loading caused 
by sodium channels amounts to about 8.5% of the bare membrane capacitance at the rest 
potential. However, it can transiently be almost 42% of the bare capacitance when one-half 
of the sodium gates is open. The membrane time constant at rest Trest = C'rcst/Grcst, where 
Grest is the corresponding membrane slope conductance, is increased accordingly. Therefore, 
the spiking dynamics slows down, cf. Fig. [T](b). 
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IV. CHANNEL NOISE AND GATING CHARGE EFFECTS 

The behavior predicted by the deterministic model should also carry implications for the 
spontaneous spiking. Especially, we expect a reduction of spiking activity due to gating 
charge effects for membrane patches of finite size. This indeed is the case as we shall 
demonstrate below. 



One of the ma 
action 



or consequences of intrinsic channel noise is the initiation of spontaneous 
quantitative measure for the occurrence of the action 



! oi tne maior consequences oi 
potentials 

mm\ oil. A 



potentials is the mean interspike interval, i.e., 

1 ^ 

iT) = j;^Y.(^^-U^i), (14) 

i=l 

where ti, i = 1, ..,N are the times for the occurrence of spike events and we set to = 0. 
Since the strength of the channel noise depends on the size of the membrane patch, (T) 
is a function of the patch size 5*. With increasing noise level or decreasing patch sizes 
S, respectively, the spike production increases and thus the mean interspike interval (T) 
decreases and even can approach the refractory time jisl, 

The spike occurrences ti are extracted from the voltage train V{t) which we obtain from 
numerical integration of the generalized HH model, cf. Eq. ffTTl) . Eq. ([5]) and Eqs. fl6al) - fl6cl) . 
The integration is carried out by a standard Euler algorithm with a step size of 1 fis. The 
Gaussian random numbers are generated by the "Numerical Recipes" routine ran2 using 



the Box-MuUer algorithm 45[ . To ensure the confinement of the gating variables between 
(all gates are closed) and 1 (all gates are open) we have implemented numerically reflecting 
boundaries at and 1. The occurrences of action potentials are determined by upward 
crossings of the membrane potential y of a certain detection threshold. Due to the very 
steep increase of membrane potential at firing, the actual choice of the detection threshold 
does not affect the results. For each trajectory we performed 2^^ simulation time steps. 
Moreover, an ensemble averaging over 100 different trajectories was done. 

In Fig. [2] we depict the dependence of the mean interspike interval (T) on the patch 
size S. Due to the inhibitory nature of the gating charge, the mean interspike interval for a 
given membrane patch is significantly increased as compared to the case where gating charge 
effects are neglected. This leads to a diminishing strength of the channel noise component, 
mostly responsible for stochastic self-excitation. To put it differently, stochastic fluctuations 
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FIG. 2: The dependence of the mean interspike interval (T) on the size S of the membrane patch is 
depicted. The accuracy of the numerical simulation results is indicated by the error bars. If gating 
currents are considered, the spontaneous spiking activity, see the grey line, is reduced compared 
to the case for which gating charge effects are neglected, indicated by the black line. Especially, 
for large patch sizes or large numbers of ion channels, gating charge effects lead to a drastically 
reduced spontaneous spiking activity. In the inset, we depict the same data on the log-log scale. 

of the membrane capacitance partially compensate the channel noise effect. Because the 
mean interspike interval is exponentially sensitive to the channel noise intensity, even a small 
relative reduction of noise intensity results in an exponentially large effect for large membrane 
patches. Therefore, the spontaneous spiking activity is drastically reduced. For example, for 
patch sizes around 45 fim^ the probability of an occurrence of a spontaneous action potential 
is one order of magnitude less in the case of gating charge effects as compared to the case 
when gating charge effects are neglected. For larger patch sizes the discrepancy becomes 
ever more striking. In clear contrast, the gating charge effects can safely be neglected for 
very small patch sizes for which the channel noise rules the dynamics. Note that even though 
in this region a discrete stochastic modeling is more appropriate, recent results show, that 
the Langevin approach does work quite well also for small patch sizes {23 1. 
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FIG. 3: Same stochastic data as in Fig. [21 but evaluated for the coefficient of variation CV, cf. 
Eq. (US}. 

V. COHERENCE OF SPONTANEOUS SPIKING ACTIVITY 



Next we investigate the influence of the gating charge on the regularity of the spontaneous 
spiking. A proper measurement is provided by the coefficient of variation, CV, which is given 
as the ratio of standard deviation to the mean value of the interspike intervals, i.e.. 



CV 



(T) 



(15) 



with the mean squared interspike interval given by (T^) := ~ ^j-i)^- For a fully 

disordered point process, meaning for a pure Poisson process, the coefficient of variation 
CV assumes the value CV = 1. For a more ordered process CV attains a smaller value, 
and for a deterministic signal it vanishes. 



Previous studies 



15 



16| where such gating current effects have been neglected, have 
demonstrated that the CV exhibits a distinct minimum for an optimal patch size S* ~ 1 fim^. 



In our case again, the spiking is mostly 



termed intrinsic coherence resonance 



m 



regu. 
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ar , cf. Fig. [31 This phenomenon has been 



47( 1 . In Fig. [3l the coefficient of variation is 



plotted vs. the patch sizes for the two cases, with and without the gating charge effect 
considered: the CV shows a distinct minimum for both cases. Though the most regular 
spiking could be observed at almost the same patch size, namely S = 1 /xm^, the spiking 
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coherence slightly deteriorates due to the gating currents' influence. The main effect of 
gating charge on stochastic dynamics is thus to slow down the spiking activity and this 
effect can be essential, cf. Fig. El 



VI. CONCLUSION AND OUTLOOK 



We extended a stochastic description of the HH model accounting for inherent channel 
noise by including gating current effects, which equivalently can be described in terms of 
capacitance fluctuations, see in eq. ( |T2l) . Our study revealed that while the deterministic HH 
model with gating charge effects does not differ dramatically from the original model for the 
standard set of parameters, the corresponding stochastic model does behave very differently 
for intermediate-to-large membrane patch sizes. A main finding is that spontaneous spiking 
activity becomes drastically reduced. The physical reason of this is that the gating current of 
translocated gating charges in sodium channels is counter-directed to the electrical current of 
sodium ions, leading finally to excitation after the whole system is driven by external current 
beyond a threshold. This results in the reduced channel noise intensity. Simultaneously, 
there occurs a rise of the membrane time constant resulting in longer spiking periods for 
supra-threshold driving. 

For the parameters studied, i.e. for the original HH parameters, the effects are relatively 
small in the noiseless case. Nevertheless, even deterministically, gating charge effects lead 



ultimate^ 
duction 



y to the existence of optimal ion channel densities for neuronal information trans- 



31 



32l | . Notably, an increase of the sodium channel density by a factor of 20 (as in 
mammalian Ranvier nodes) would cause gating charge effects of increased importance. The 
striking feature revealed in this work is that the gating charge effects can play a fundamental 
role in the spontaneous spiking dynamics of intermediate-to-large membrane patches even 
if they are still relatively minor deterministically. The reduction of spontaneous spiking 
exceeds one order of magnitude starting from a patch size around 45 fim^ for the parameters 
studied. This provides an inherent mechanism for the channel noise reduction in neurons. 
This might explain the discrepancy between theoretical predictions of stochastic HH model 
and some experimental results for real neurons, e.g. see Ref. |48|. However, the effect 
of intrinsic coherence resonance as signature of the influence of channel noise in excitable 
membranes still remains. 
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We share the confident behef that our investigations of the influence of gating charge 
effects on the channel noise-induced spiking activity in an archetypal model provide some 
new insights into the underlying physical principles and mechanism of neuronal signaling. 
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Glossary 

Coherence Resonance. Noise-induced improvement of the regularity of the system's output. 
If the noise stems from internal dynamics the effect is called intrinsic coherence reso- 
nance. 

Hopf-bifurcation. A change from a topology with flxed point solution to a topology with 
oscillatory solution under a small variation of a parameter within a nonlinear system. 

Ito-Stratonovich dilemma. Interpretation problem which arises in the context of Langevin 
equations in case of multiplicative Gaussian white noise. 

Kramers-Moyal expansion. Procedure to transform a master equation for a discrete stochas- 
tic process into a Fokker-Planck equation for approximate continuous stochastic pro- 
cess. 

Poisson process. A stochastic process which is memory-less (Markovian) with exponentially 
distributed waiting times between two successive events. 

Stochastic Resonance. A anomalous, noise-assisted enhancement of transduction of weak 
(deterministic or stochastic) signals. 
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